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ABSTRACT 



<—^ Aims. We aim to study the physical conditions, particularly the excitation state, of the intermediate-temperature gas components in the protoplan- 
—5 , etary nebula CRL 618. These components are particularly important for understanding the evolution of the nebula. 

Methods. We performed Herschel/HIFI observations of several CO lines in the far-infrared/sub-mm in the protoplanetary nebula CRL 618. The 
high spectral resolution provided by HIFI allows measurement of the line profiles. Since the dynamics and structure of the nebula is well known 
from mm-wave interferometric maps, it is possible to identify the contributions of the different nebular components (fast bipolar outflows, double 
shells, compact slow shell) to the line profiles. The observation of these relatively high-energy transitions allows an accurate study of the excitation 
conditions in these components, particularly in the warm ones, which cannot be properly studied from the low-energy lines. 
Results. The 12 CO 7=16-15, 10-9, and 6-5 lines are easily detected in this source. 13 CO 7=10-9 and 6-5 are also detected. Wide profiles 

C"j showing spectacular line wings have been found, particularly in 12 CO 7=16-15. Other lines observed simultaneously with CO are also shown. 

Q,, Our analysis of the CO high-7 transitions, when compared with the existing models, confirms the very low expansion velocity of the central, 
I ■ dense component, which probably indicates that the shells ejected during the last AGB phases were driven by radiation pressure under a regime 

O ' of maximum transfer of momentum. No contribution of the diffuse halo found from mm-wave data is identified in our spectra, because of its low 
temperature. We find that the fast bipolar outflow is quite hot, much hotter than previously estimated; for instance, gas flowing at 100 km s must 



OS 



5-H 



CZ3 . have a temperature higher than ~ 200 K. Probably, this very fast outflow, with a kinematic age < 100 yr, has been accelerated by a shock and has 
■ not yet cooled down. The double empty shell found from mm-wave mapping must also be relatively hot, in agreement with the previous estimate. 

Key words, stars: AGB and post-AGB - stars: circumstellar matter, mass-loss - planetary nebulae - planetary nebulae: individual: CRL 618 

>: 

[^. \ 1 . Introduction in molecular lines dBujarrabal et al.ll200lT) . PPNe have been ac- 

\f\ . curately observed in mm-wave lines, particularly by means of 

^ , Protoplanetary nebulae (PPNe) are known to present very fast inter ferometric maps with resolutions ~ 1". Thanks to those ob- 

bipolar outflows, along with slower components, which are serV ations, the structure, dynamics, and physical conditions in 

I probably the remnants of the mass ejection during the previous these nebulae ^ often quite well known HoweV er, observa- 

^ , AGB phase. The bipolar flows typically reach velocities of 100 tions of me low _ y transitions are not very use ful for studying 

^ ; km s \ and affect a sizablefraction of the nebular mass, ~ 0.1 - the warm gas componen ts. The well-studied 7=2-1 and 7=1-0 

. . . 0.3 M ( |Buiarrabaletal. || 2001|) . These dense flows actually rep- transitions of i2 CO only require tem peratures of T k ~ 15 K to 

m > ■ resent intermediate states in the spectacular evolution from the be exdted Inde ed their maximum emissivity occurs for excita- 

kj ; spherical and slowly expanding circumstellar envelopes around tion tem p er atures of 10 - 20 K, and the line intensities and line 

^ , AGB stars to the planetary nebulae, which usually show bipo- intens i ty ratios on i y slightly depend on the excitation state in rel- 

d ■ lar or ring-like symmetries. Such remarkable dynamics is proba- atively warm gas Needless t0 say? observations in the visible or 

bly the result of the interaction between the AGB and post-AGB near infrared ranges tend to select hot regions> with typical tem . 

winds: axial, very fast post-AGB jets colliding with the denser per atures over 1000 K. The proper study of warm regions, 100 

material driven isotropically away from the star during its AGB K < ^ < 1000 K therefore re q U i re s observations at intermediate 

phase (e.g. [ Bahck & Frank| l 2002|) . The presently observed bipo- wave i e ngths, in the far infrared (FIR) and sub-mm ranges. 

lar outflows would then correspond to a part of the relatively 

dense shells ejected during the last AGB phase, mostly their po- Because of the role of shocks in PPN evolution, these 
lar regions, accelerated by the shocks that propagate during the warm regions are particularly important for understanding neb- 
PPN phase. ular structure and e volution. In some well-studied cases, e.g. 
The massive bipolar outflows in PPNe, as well as the unal- M 1-92 and M 2-56 (Alcol ea et al.ll2007llBuiarrabal et all 19981 
tered remnants of the AGB shells, usually show strong emission ICastro-Carrizo et al.ll2002l) . the high-velocity, massive outflows 
are found to be very cold, with temperatures < 20 K, which im- 



* Herschel is an ESA space observatory with science instruments pro- phes very fast cooling in the shock-accelerated gas. No warm 
vided by European-led Principal Investigator consortia and with impor- component representing the gas recently accelerated by the 
tant participation from NASA. shock front has been identified in these sources. In other cases 
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Fig. 1. HIFI observations containing detected 12 CO and 13 CO 
lines in CRL 618 {T m \, vs. rest frequencies). Frequency scales 
and detected lines are indicated; note that the observations are 
performed in DSB mode. 



such as CRL 618, interferometric imaging of the 12 CO 7=2-1 
line show s that dense gas in axial st ructures presents higher tem- 
peratures (San chez Contreras etai1l2004l hereafter, SC04). But, 
precisely because of their relatively high excitation, the temper- 
ature estimate in these components from CO 7=2-1 is very un- 
certain. Indeed, even the presence of such high excitation in the 
dense bipolar outflows in CRL 618 remained to be demonstrated. 

Other atte mpts to study the warm gas in CRL 618 were car- 
ried out by (i) lJusttanont et a l. (2000) from ISO data with low 
spectral resolution; (;';") IPardo et alJ (120041) . w ho focused on the 
chemi stry of the different components; and (hi) Nakashim a et alJ 
(120071) . who also obtained maps of the 7=6-5 transition in 
CRL 6 18, but with less detail than in SC04. 

The Herschel Space Telescope is well-suited to studying 
warm gas around evolved stars in the FIR and sub-mm. The high 
spectral resolution that can be achieved with its heterodyne in- 
strument HIFI (better than lkms -1 ) is particularly useful for 
this purpose, since kinematics offers a fundamental key to under- 
standing this warm, shocked gas. Here we present Herschel/HIFI 
observations of CRL 618 in several molecular lines of I2 CO and 
13 CO that were obtained as part of the guaranteed-time key pro- 
gram HIFISTARS, which is devoted to the study of intermediate- 
excitation molecular lines in nebulae around evolved stars. 



Table 1. Summary of the Herschel telescope characteristics, 
line band DSB r sys HPBW T mb /T. d cal. uncert. 
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2. Observations 



We used the Hersch el/HIFI instrument dPilbratt et alJ 1201 Ol 
de Graauw et al. 201.3) to observe the 7=6-5, 10-9, and 16-15 
transitions of 12 CO and 13 CO in the PPN CRL 618 ( 13 CO 
7=16-15 was not detected); see Figs. 1 and 2. Other molecular 
lines were also detected within the observed frequency ranges. 
The data were taken using the two orthogonal HIFI receivers 
available at each band. Both receivers work in double side-band 
(DSB) mode, which effectively doubles the instantaneous IF 
coverage. Care was taken when choosing the local oscillator fre- 
quency, maximizing the number of observed interesting lines. 

The observations were obtained in the dual-beam-switching 
(DBS) mode. In this mode, the HIFI internal steering mirror 
chops between the source position and an emission-free position 
3' away. The telescope then alternatively locates the source in 
either of the chopped beams, providing a double-difference cal- 
ibration scheme, which allows a more efficient cancellation of 
the residual standing waves in the spectra (see additional details 
in Helmich et al. 2010). This procedure works very well except 
for the 7=16-15 lines, where strong ripples were found in some 
spectra, especially in the V-receiver. 

The HIFI data shown here were taken using the Wide-Band 
Spectrometer (WBS), an acousto-optical spectrometer that pro- 
vides simultaneous coverage of the full instantaneous IF band 
in the two available orthogonal receivers, with a spectral reso- 
lution of 1 . 1 MHz. The data shown in all the figures have been 
resampled and smoothed to a resolution of about 2 km s _1 . 

The data were processed with the standard HIFI pipeline 
using HIPE, with a modified version of the level 2 algorithm 
that yields unaveraged spectra with all spectrometer sub-bands 
stitched together. Later on, the spectra were exported to CLASS 
using the hiClass tool within HIPE, for further inspection, flag- 
ging out data with outstanding ripple residuals, final averaging, 
and baseline removal. We checked that the profile of the 12 CO 
7=16-15 line, in particular, is not significantly affected by rip- 
ples, which in fact are not noticeable even in the relatively flat 
parts of the final spectrum. The data were originally calibrated 
in antenna temperature units and later converted into main-beam 
temperatures (7 m b). In all cases we assumed a side-band gain 
ratio of one. A summary of the telescope characteristics and ob- 
servational uncertainties is given in Table 1 . 

We also report here 12 CO 7=4-3 and 7-6 lines observed 
from the ground with the APEX telescope. These observa- 
tions were performed in preparation for the HIFI observations 
in 2006 and more details will be given in a future paper. We 
used the FLASH receiver equipped with two FTS spectrometers, 
which allows simultane ous observation of the two CO lines; see 
Hevminc k et alJ (120061) for a description of the system. The re- 
sulting APEX spectra are shown in Fig. 2, after being rescaled 
to 7mh units using the values for the telescope efficiencies given 
bv lGusten et al.l (12006b . In this figure, we also show IRAM 30m 
data data for 12 CO 7= 1 -0 and 2- 1 from lBuiarrabal et al.l(l200lh . 
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3. Nebula model 

Detailed mapping of CO emission at 1mm wavelength was per- 
formed by SC04, who derived the physical conditions, structure, 
and dynamics of the nebula from model fitting of their maps. 
Several warm components were identified (see Fig. 3): a com- 
pact dense core with temperatures ~ 100 K (decreasing with dis- 
tance to the center), a double (empty) shell with a typical temper- 
ature of 200 K, and a very fast bipolar outflow (running inside 
the cavities) with temperatures 5 100 K (also decreasing with 
distance). Another cooler component was found: a diffuser halo 
expanding at 17 km s -1 . The empty shells show a roughly el- 
liptical shape that strongly suggests that they are the result of 
a bow-like shock, but with a moderate expansion velocity. The 
expansion velocity of the compact component was found to be 
particularly low, 5 10 km s _1 . 

As a starting point, we adopted a similar description of the 
nebula structure. Calculation of the expected emission in our 
case must be more sophisticated, because the high-/ CO transi- 
tions cannot be assumed to be thermalized in the whole nebula. 
Therefore, we determined the CO level population performing 
LVG calculations for a large number of points in the model neb- 
ula. LVG calculations are fully justified here because of the high 
velocity gradients characteristic of this source. With these values 
of the populations we calculated the emissivity and absorption 
coefficients at each point. Finally, the brightness distribution is 
calculated along a number of lines of sight and is convolved with 
the telescope beam shape, described by a Gaussian function. The 
results are profiles in units of main-beam temperature, directly 
comparable to our observations. Therefore, opacity effects are 
properly taken into account, both in the excitation and line pro- 
file calculations. We have checked that the high-/ lines are often 
underexcited in our case, particularly for densities under ~ 10 6 
cirT 3 . 

More detailed discussions of the code and involved uncer- 
tainties will be presented in a forthcoming paper. Here we mostly 
discuss a preliminary comparison of the model results with our 
observations of the I2 CO /=16-15 line. The model fitting can- 
not be considered as satisfactory without studying all the ob- 
served profiles, as well as the mm-wave maps, requiring a very 
careful treatment of the numerical uncertainties and a detailed 
discussion of the nebula structure. 



4. Results 

We present in Fig. 2 our observational results for the 12 CO lines 
obtained with HIFI, together with other lines observed from the 
ground (Sect. 2). The baseline of the /=16-15 line is not well 
determined, because of the lack of spectral coverage and the 
presence of other lines, which mostly affects the intensity of the 
line wings at extreme velocities. 

We have mentioned (Sect. 1) the interest of observing high-/ 
transitions in order to better estimate the excitation conditions in 
warm regions. Our highest transitions require several hundred K 
to be significantly excited; for instance, the energy of the 12 CO 
J =16 level is equivalent to 750 K. 

We have compared our data with the predictions of the model 
presented in Sect. 3, originally developed to explain the mm- 
wave maps (see assumed nebula structure in Fig. 3). This model 
can reasonably explain the central part of the high-/ emission 
with only moderate changes. We can see in Fig. 4 the compar- 
ison between the observations and predictions for a model in 
which the density and temperature of the shells have been in- 
creased by 25% (dashed, red line). The asymmetry in the ob- 



CO J=1 6—1 5_ 

(HIFI) 




Fig. 2. 12 CO lines observed in CRL 618 with HIFI, in units of 
r m b vs. LSR velocity. Other CO lines observed from the ground 
are also shown. Note that the high-velocity wings of the 12 CO 
/=10-9 profile are affected by emission of other lines (Fig. 1). 



served profile cannot be explained by radiative transfer phenom- 
ena, and it probably reveals true asymmetries with respect to the 
equator that are not included in our models. The central spectral 
feature is very narrow and comes from the low-velocity compact 
component of the nebula; in fact, the fitting is improved if the 
velocity of central component is still decreased, to typical values 
< 5 km s , but keeping the strong velocity gradient characteris- 
tic of this region, in order to reproduce the triangular shape of 
the observed feature. The two intense humps at both sides of the 
central maximum would come from the hollow shells. 

However, the high-velocity wings are severely underesti- 
mated by our standard model, and we would have to significantly 
increase the excitation conditions of the very fast bipolar outflow 
of CRL 618 to reproduce their intensity. We can also see in Fig. 
4 our predictions for a model similar to the previous one, but 
with 7\ ~ 200 K in the regions that present expansion velocities 
~ 100 km s -1 . Other parameters of the fast outflow, particularly 
its velocity and density distributions do not change with respect 
to the original model. The asymmetry between the red and blue 
line wings is not reproduced by our predictions also in this case. 
It is remarkable that the high-velocity outflow obviously con- 
tributes to the emission at the profile central features. Therefore, 
in this model the requirements to reproduce the secondary max- 
ima are significantly weaker, and a temperature similar to or even 
lower than assumed in our original model for the empty shells 
would be compatible with the observations. The general proper- 
ties of the central, dense component mentioned before, in partic- 
ular its low velocity, are required in this case too. The high ve- 
locity wings are also detected in emission from other molecules, 
such as H2O, HCN, and CN (Fig. 1), which must be significantly 
abundant in this recently shocked gas. 
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Fig. 3. Basic model properties used in our calculations, see 
Sects. 3 and 4 (adapted from SC04). 




CO J=I6-15 



Fig.4. 12 CO 7=16-15 observed line and predictions of our 
model (Sects. 3, 4), without (dashed, red line), and with (dot- 
dashed, green line) significant increase of the excitation in the 
very fast outflow. 



The relatively high temperatures deduced here for the fast 
bipolar flow relax the discrepancy usually found between the 
high excitation of the shocked gas predicted by theoretical mod- 
els and the observational results (an i ntricate theoreti cal problem 
not discussed in this lett er, see e.g . Lee et alj|2009l) . However, 
the discrepancy persists. iLee et all (120091) predict temperatures 



of the high-velocity gas in CRL 618 over ~ 1000 K and too weak 
CO emission in all rotational lines, since shocks are expected to 
dissociate molecules. 

The low-excitation component of the nebula model by SC04, 
the extended halo, has apparently no counterpart in the observa- 
tions, because the model predicts a very low intensity from such 
cool gas and the high-./ line profiles do not seem to require any 
contribution from it. 

In Fig. 1 we also show our HIFI spectra of 13 CO lines. The 
13 CO 7=16-15 observations are not very sensitive, so the line is 
not detected with a limit T m \, < 0.2 K. As we can see, the contrast 
between 12 CO and 13 CO lines is high (mainly in the line wings, 
about a factor ten for the highest transitions). This result is com- 
patible with our calculations, which suggest moderate opacities 
in high-7 12 CO lines from the main nebular components, in par- 
ticular with t(16-15) < 1 for gas flowing at > 100 km s . 

Our results can therefore be summarized as follows 



1. We detected high-7 CO emission using Herschel/HIFI. The 
high-velocity line wings characteristic of this source become 
progressively dominant as the level energies increase, with 12 CO 
7=16-15 showing a spectacular composite profile. 

2. The temperature of the very fast bipolar outflow in CRL 618 
is high, significantly higher than the previously adopted val- 
ues. SC04 proposed a temperature for this component < 100 K, 
which, in view of the intense line wings seen in the 7=16-15 
transition, must be significantly increased. From our calcula- 
tions, we estimate that gas flowing at about 100 km s -1 must 
have a temperature ~ 200 K. We suggest that this very fast out- 
flow, with a kinematic age < 100 yr, was accelerated by a shock 
and has not yet fully cooled down. The rest of the physical con- 
ditions are not significantly changed, therefore the dynamical 
parameters (including the high momentum and kinetic energy) 
remain the same as those deduced by SC04. 

3. We confirm the low expansion velocity of the very dense, cen- 
tral component. This low velocity coincides with a significant 
increase in the mass-loss rates during the last ~ 400 yr of the 
AGB phase (SC04). We suggest that the significant decrease in 
the expansion velocity is caused by the ejection of material by 
the star during the last AGB phases being driven by radiation 
pressure under a regime of maximum momentum transfer from 
radiation to gas. 

We recall that our analysis is preliminary. Further develop- 
ments, in particular including fitting the mm-wave maps and our 
HIFI data, are required. 

Acknowledgements. HIFI has been designed and built by a consortium of 
institutes and university departments from across Europe, Canada, and the 
United States under the leadership of SRON Netherlands Institute for Space 
Research, Groningen, The Netherlands, and with major contributions from 
Germany, France, and the US. Consortium members are: Canada: CSA, 
U.Waterloo; France: CESR, LAB, LERMA, IRAM; Germany: KOSMA, 
MPIfR, MPS; Ireland, NUI Maynooth; Italy: ASI, IFSI-INAF, Osservatorio 
Astrofisico di Arcetri- INAF; Netherlands: SRON, TUD; Poland: CAMK, CBK; 
Spain: Observatorio Astronomico Nacional (IGN), Centro de Astrobiologi'a 
(CSIC-INTA). Sweden: Chalmers University of Technology - MC2, RSS & 
GARD; Onsala Space Observatory; Swedish National Space Board, Stockholm 
University - Stockholm Observatory; Switzerland: ETH Zurich, FHNW; USA: 
Caltech, JPL, NHSC. HCSS / HSpot / HIPE is a joint development (are joint de- 
velopments) by the Herschel Science Ground Segment Consortium, consisting 
of ESA, the NASA Herschel Science Center, and the HIFI, PACS, and SPIRE 
consortia. This work has been partially supported by the Spanish MICINN, pro- 
gram CONSOLIDER INGENIO 2010, grant "ASTROMOL" (CSD2009-00038). 
R.Sz. and M.Sch. acknowledge support from grant N 203 393334 from the 
Polish MNiSW. K.I. acknowledges the funding from SNSB. J.C. acknowledges 
funding from MICINN, grant AYA2009-07304. This research was performed, 
in part, through a IPL contract funded by the National Aeronautics and Space 
Administration. 



References 

Alcolea, J., Neri, R., & Bujarrabal, V. 2007, A&A, 468, L41 

Balick, B., & Frank, A. 2002, ARA&A, 40, 439 

Bujarrabal, V., Alcolea, J., & Planesas, P. 1992, A&A, 257, 701 

Bujarrabal, V., Alcolea, J., & Neri, R. 1998, ApJ, 504, 915 

Bujarrabal, V., Castro-Carrizo, A., Alcolea, J., & Sanchez Contreras, C. 2001, 

A&A, 377, 868 
Castro-Carrizo, A., Bujarrabal, V., Sanchez Contreras, C, Alcolea, J., & Neri, 

R. 2002, A&A, 386, 633 
de Graauw, Th., et al. 2010, A&A, in press 
Giisten, R., Nyman, L. A., Schilke, P., Menten, K., Cesarsky, C, & Booth, R. 

2006, A&A, 454, L13 
Helmich, E, et al. 2010, A&A, this volume 
Heyminck, S., Kasemann, C, Glisten, R., de Lange, G., & Graf, U. U. 2006, 

A&A, 454, L21 
Justtanont, K., et al. 2000, A&A, 360, 1117 
Lee, C.-E, Hsu, M.-C, & Sahai, R. 2009, Apl, 696, 1630 
Nakashima, J.-i., et al. 2007, Al, 134, 2035 
Pardo, I. R., Cernicharo, I., Goicoechea, J. R., & Phillips, T. G. 2004, ApJ, 615, 

495 



V. Bujarrabal et al.: Herschel/HIFI observations of high-./ CO transitions in the protoplanetary nebula CRL 618 

Pilbratt, P., et al. 2010, A&A, in press 

Sanchez Contreras, C, Bujarrabal, V., Castro-Carrizo, A., Alcolea, J., & Sargent, 
A. 2004, ApJ, 617, 1142 (SC04) 



1 Observatorio Astronomico Nacional (IGN), Ap 112, E-28803 
Alcala de Henares, Spain 

e-mail: v.bujarrabal@oan.es 

2 Observatorio Astronomico Nacional (IGN), Alfonso XII N°3, E- 
28014 Madrid, Spain 

3 CAB, INTA-CSIC, Ctra de Torrejon a Ajalvir, km 4, E-28850 
Torrejon de Ardoz, Madrid, Spain 

4 Instituut voor Sterrenkunde, Katholieke Universiteit Leuven, 
Celestijnenlaan 200D, 3001 Leuven, Belgium 

5 Sterrenkundig Instituut Anton Pannekoek, University of 
Amsterdam, Science Park 904, NL-1098 Amsterdam, The 
Netherlands 

6 Onsala Space Obseravtory, Dept. of Radio and Space Science, 
Chalmers University of Technology, SE^13992 Onsala, Sweden 

7 European Space Astronomy Centre, ESA, P.O. Box 78, E-28691 
Villanueva de la Canada, Madrid, Spain 

8 Harvard-Smithsonian Center for Astrophysics, Cambridge, MA 
02138, USA 

9 Max-Planck-Institut fur Radioastronomie, Auf dem Hiigel 69, D- 
53121 Bonn, Germany 

10 The Johns Hopkins University, 3400 North Charles St, Baltimore, 
MD 21218, USA 

11 Department of Astronomy, AlbaNova University Center, 
Stockholm University, SE-10691 Stockholm, Sweden 

12 Joint ALMA Observatory, El Golf 40, Las Condes, Santiago, 
Chile 

13 N. Copernicus Astronomical Center, Rabiahska 8, 87-100 Toruh, 
Poland 

14 Department of Astrophysics/IMAPP, Radboud University 
Nijmegen, Nijmegen, The Netherlands 

15 Astronomical Institute, Utrecht University, Princetonplein 5, 3584 
CC Utrecht, The Netherlands 

16 Instituto de Radioastronomia Milimetrica (IRAM), Avda. Divina 
Pastora 7, E-18012 Granada, Spain 

17 Observatorio Astronomico Nacional (IGN), Centro Astronomico 
de Yebes, Apartado 148, E-19080 Guadalajara, Spain 

18 SRON, Netherlands Institute for Space Research, Landleven 12, 
9747 AD Groningen, The Netherlands 

19 Sterrewacht Leiden, University of Leiden, P.O. Box 9513, 2300 
RA Leiden, The Netherlands 

20 KOSMA, I. Physik. Institut, Universitat zu Koln, Zulpicher Str. 
77, D 50937 Koln, Germany 



